Formation of atom wires on vicinal silicon 
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The formation of atomic wires via pseudomorphic step-edge decoration on vicinal silicon surfaces 
has been analyzed for Ga on the Si(112) surface using Scanning Tunneling Microscopy and Density 
Functional Theory calculations. Based on a chemical potential analysis involving more than thirty 
candidate structures and considering various fabrication procedures, it is concluded that pseudo- 
morphic growth on stepped Si(112), both under equilibrium and non- equilibrium conditions, must 
favor formation of Ga zig-zag chains rather than linear atom chains. The surface is non-metallic 
and presents quasi-one dimensional character in the lowest conduction band. 

PACS numbers: 68.35.-p, 68.37.-d, 81.07.-b, 73.20.-r 
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Deposition of submonolayer amounts of metal ad- 
sorbates on high-index silicon surfaces often results in 
the formation or "self assembly" of single-domain chain 
structures or "atomic wire arrays" 0. The basic idea 
is that metal atoms should preferentially adsorb at the 
highly reactive step-edges and form parallel chains of 
metal atoms. The transport properties of atomic chains 
are currently of great interest because the venerable 
quasi-particle concept of Fermi liquid theory is expected 
to break down in one dimension. Photoemission from 
atomic gold chains on Si(557) |2| and Si(553) [3| re- 
vealed dispersive metallic bands that are strongly one- 
dimensional but these surfaces turned out to be compli- 
cated reconstructions with large unit cells. The metal 
atoms are being incorporated into the terraces rather 
than at the step-edges and the observed metallicity is be- 
lieved to originate from the substrate's dangling bonds. 
Surface reconstructions generally represent a quirk of na- 
ture and it therefore remains difficult to predict or engi- 
neer the formation of useful structures. 

A precise definition of interfaces at the atomic scale is 
usually best achieved by epitaxial or pseudomorphic (i.e. 
strained overlayer) growth. In this Letter, we investigate 
the general feasibility of atom wire formation by pseu- 
domorphic step-edge decoration for the case of Si(112) 
with Ga adatoms 0, 0, 0, • As with any surface recon- 
struction, the geometrical configuration of the atom wires 
should reflect the balance between energy gain through 
passivation of dangling-bond surface states, and energy 
cost from the resultant increase of bond angle strains. In 
the case of pseudomorphic overlayers, one rather speaks 
of misfit strain, which can be relieved by forming "misfit 
dislocations" or, in the case of a monoatomic overlayer, 
vacancy-line superstructures. 

The general question that has not been addressed is 
whether atomic chain structures can actually represent 



stable or metastable structures or whether additional 
metal atoms should passivate the dangling bonds of the 
terraces, thereby forming a more complex overlayer. We 
present new Scanning Tunneling Microscopy (STM) data 
of Si(112)6 x 1-Ga that reveal unprecedented detail, in- 
cluding a hitherto unnoticed symmetry breaking in the 
vacancy lines, ruling out the formation of the linear atom 
chains that had been proposed before 0, 0, S @- 
The STM data are complemented with extensive Den- 
sity Functional Theory (DFT) calculations and theoreti- 
cal STM simulations. Our analysis shows that Ga atoms 
adsorb on the terraces as well as at the step-edges, form- 
ing double chains or zig-zag chains. A chemical potential 
analysis confirms that the zig-zag structure represents 
the thermodynamically stable structure. The passiva- 
tion of the substrate's dangling bonds clearly dictates 
the structure of the metal overlayer; the alternative of 
creating buckling distortions on the terraces to eliminate 
partially-filled dangling bonds appears less effective 
in lowering the total energy, which precludes the forma- 
tion of single-atom wires. This simple passivation ar- 
gument should have general validity for semiconductor 
interfaces, which suggests that single-atom linear wires 
cannot be stabilized under the quasi-equilibrium condi- 
tions of a molecular beam epitaxy experiment, and other 
alternatives (e.g. zig-zag chains) must be explored. 

Experiments were carried out in an ultra-high vacuum 
system (base pressure < 5 x 1Q — 11 mbar) with a Ga 
deposition source, direct current sample heating facili- 
ties and a variable temperature STM. The sample was 
flashed at 1475 K to remove the native oxide. During 
resistive heating, the current was directed parallel to the 
nano-facets of the clean (112) surface (i.e. in the [110] 
direction) in order to avoid current-induced step bunch- 
ing. Ga was deposited using a commercial effusion cell. 
We have prepared the Si(112)6 x 1-Ga surface in two 
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different ways. In the "one-step" procedure, Ga is evap- 
orated onto a Sifl 12) substrate kept at a temperature 
of T s — 825 K |. In the "two-step" procedure, Ga is 
deposited onto a Si(l 12) substrate kept at room temper- 
ature. After depositing more than 1 ML, the substrate is 
annealed to T s = 825 K to remove excess Ga and form the 
6x1 structure . Experimentally, the one-step and two- 
step growth procedures always produced identical 6x1 
structures. STM experiments were performed at room 
and at low (~ 40 K) temperature using etched tungsten 
tips. STM images of the filled and empty electronic states 
were obtained with a constant current between 0.06 and 
0.2 nA and bias voltages between 1 and 2 V. 




FIG. 1: (a) Empty state STM image of the Si(112)6 x 1-Ga 
surface. Tunneling conditions: 1.5 V, 0.2 nA. (b) Filled state 
STM image. Tunneling conditions: -2.0 V, 0.1 nA. (c) Empty 
state and (d) filled state theoretical STM images correspond- 
ing to surface structure of Fig. [3] 

Fig.^a) shows an atomic resolution empty state STM 
image of the Si(112)6 x 1-Ga surface. Two parallel atomic 
lines are observed per unit cell, running in the [110] direc- 
tion, intersected by quasi-periodical vacancy lines. The 
brightest atomic lines in Fig. [Ha) are the same atomic 
lines as imaged by Baski et al. |6( . But in addition we ob- 
serve an atomic line lying in between the brighter lines. 
These two parallel atomic lines appear to form a zig- 
zag pattern, resulting in a zig-zag asymmetry in the va- 
cancy line. In the filled state image, Fig. QI&), a rela- 
tively bright profusion is the dominant feature. Registry 



aligned dual bias images show that this bright protru- 
sion is located in the vacancy lines. Bright protrusions 
in different unit cells are connected by fuzzy lines along 
the [110] direction. These STM measurements show that 
the step-edge decorated structural model U, la, |y, LD, Ig 01 
the Si(112)6 x 1-Ga surface cannot be correct since the 
vacancy lines in that model have mirror plane symme- 
try (i.e. mirror plane perpendicular to [110]). We fur- 
thermore notice that low-temperature STM (40 K) does 
not provide evidence for a buckling distortion of the Si 
terrace atoms, as predicted for the step-edge decorated 
structure @. Rutherford Backscattering Spectrometry 
experiments indicate 9 ± 1 Ga atoms per 6x1 unit cell, 
as compared to 5 Ga atoms for the step-edge decorated 
model. 
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FIG. 2: Surface energy as a function of fiGa for the most 
promising Si(112)6 x 1-Ga candidate structures. The struc- 
ture with 10 Ga atoms per 6x1 unit-cell (see Fig. |2J is used 
as reference. Inset: analysis of the stability of the zig-zag 
structure (Fig. with respect to the longitudinal periodic- 
ity. Note the different vertical scale (also in eV). 

Extensive DFT calculations were performed in order to 
identify the precise atomic structure of the Ga/Si(112) 
surface. New candidate structural models correspond- 
ing to Ga-coverages ranging from 5 to 11 Ga atoms per 
6x1 unit cell were explored using an efficient local- 
orbital DFT method 0- In total, more than 30 new 
structures were fully relaxed, and their relative stabil- 
ity was analyzed comparing their surface energies F — 
Etot — ^SiNsi — ^GaNca ■ since different surface struc- 
tures contain different numbers of Ga and/or Si atoms 
we have to consider also the Ga and Si chemical poten- 
tials noa and usi {Etot is the total energy per 6x1 
unit-cell, and N$i, Nca ar e the number of Si and Ga 
atoms per 6x1 unit-cell). For fisi we have used the to- 
tal energy of bulk-Si, and we have analyzed the stability 
of the different surface structures as a function of UGa- 
The most promising surface structures were further re- 
fined using a plane- wave (PW) DFT code Figure 
121 shows the surface energy as a function of [iQa-, for the 
most stable surface structures containing 5, 6, 9, 10 and 
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11 Ga atoms ^lj- The chemical potential of bulk-Ga, 
ua a {bulk), can be expected to be the upper limit for fiQa 
[6j . The actual chemical potential [iQ a of the overlayer 
can be estimated from the experimental conditions for 
Ga deposition. Consider the one-step process in which 
the (6 x l)-phase is formed under a Ga flux at a sample 
temperature of T s = 825 K. At this temperature the in- 
coming flux is balanced by a flux of Ga atoms desorbing 
from the 6x1 surface and the sample can be considered 
to be in thermodynamic equilibrium. Approximating the 
chemical potential in the effusion cell by the total energy 
of bulk-Ga, p,G a {bulk), the overlayer chemical potential 
can be determined from 

MGa = ^Ga{bulk) - k B Tln(—) 

Ps 

where p c is the Ga vapor pressure in the effusion cell and 
p s the Ga vapor pressure at the sample. We obtain ^1 
VGa = PGa-bulk ~ 0.32(0.48) eV. 
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FIG. 3: Ball and stick representation of the energy minimized 
structure (zig-zag model) for the Si(112)6 x 1-Ga surface. Top 
view and side view. 

In Figure |3 we observe that the structure with 10 Ga 
atoms per 6x1 unit-cell is the most stable for this range 
of UGa- As shown in Fig. [31 the Ga atoms in this surface 
form two parallel rows: a row of 'step-edge'-Ga atoms 
adsorbed at the (lll)-like step and a second row of Ga 
atoms adsorbed at the terraces. The two Ga rows form 
a zig-zag pattern. Each Ga-row contains 5 Ga atoms per 



6x1 unit-cell, i.e. there is a Ga- vacancy in each row. In- 
side the vacancy line a rebonding has taken place, leading 
to the formation of a Si-Ga dimer at the step-edge and 
a Si-Si dimer on the terraces. In this zig-zag structure 
all the partially-filled Si dangling bonds are replaced by 
empty Ga dangling bonds, and the surface is semicon- 
ducting. Figure 0] shows the calculated band-structure 
for this surface: the Ga atoms induce a surface band in 
the upper part of the band gap, with the minimum at 
the X-point. The dispersion of this band close to the 
minimum yields an effective mass of m* ~ 1.48 along 
the X — r direction, while m* ~ 0.15 along the X — K 
direction, indicating quasi-one-dimensional conduction. 
This Ga-band is initially empty but could be populated 
on heavily-doped n-type subtrates or perhaps by appli- 
cation of an external bias. 

The step-edge decorated Si(112)6 x 1-Ga structure 
might be stabilized for very low nca values (5 Ga case in 
Fig. In order to analyze this possibility we have calcu- 
lated the total energy of a clean Si(112) surface 0] and 
compared the step-edge decorated Si(112)6 x 1-Ga case 
with a hypothetical half-and-half surface: half of the sur- 
face with the zig-zag structure ( i.e. both step-edge and 
terrace Ga rows) and the other half with clean Si(112). 
These two cases present the same average Ga coverage. 
We find the half-and-half case to be lower in energy by 
~ 0.7 eV/(10 Ga atoms) for all nca, showing that the 
step-edge decorated case is unlikely to be stabilized at 
lower coverage. We also investigated the possibility that 
the metastable step-edge decorated structure is stabilized 
kinetically through selective desorption of the Ga terrace 
atoms during the annealing step of the two-step process 
(MGa -C HGa{bulk)). Total energy calculations, however, 
show that the desorption energies for the Ga atoms of 
Fig. are lower for step-edge atoms than for terrace 
atoms by ~ 0.7 eV/atom. Hence, there seems to be no 
kinetic pathway to achieve the step-edge decorated struc- 
ture via thermal processing. 

The stability of the zig-zag model (Fig. EJ) as a function 
of the longitudinal periodicity, i. e. changing the number 
of Ga atoms between vacancy lines in each Ga row from 4 
(5 x 1) to 6 (7 x 1), is analyzed in the inset of Fig. [21 For 
our range of fj,Qa the 6x1 periodicity is the most stable, 
while the 5 x 1 is only 0.1-0.2 eV/(6 x 1 unit-cell) higher, 
in good agreement with the experimental observation of 
coexisting 5x1- and 6 x 1-like unit cells on the surface 
The preferred 6x1 periodicity of the vacancy lines 
appears identical to that of the simple Frcnkcl-Kontorova 
prediction for one-dimensional chains Q. 

Figs. W[c-d) show simulated STM images for the 
zig-zag model with a top view of the structural model 
superimposed. The agreement between the experimental 
and theoretical STM images is excellent. It shows that 
the two atomic lines imaged in the empty state STM im- 
ages are indeed the step-edge Ga line and the terrace 
Ga line. Furthermore, the zig-zag asymmetry in the va- 
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cancy lines observed experimentally is neatly reproduced 
by the simulated STM images. In the filled state image, 
fuzzy lines with a bright protrusion inside the vacancy 
line are observed, also in agreement with the experimen- 
tal images. The bright protrusion originates from the 
Si-Ga dimer in the vacancy lines, while the fuzzy lines 
are due to Si-Ga bonds on the (lll)-like terrace, close to 
the step-edge. 



FIG. 4: Surface band structure of Si(112)6 x 1-Ga around the 
semiconductor energy gap. The quasi-one dimensional state 
discussed in the text is highlighted. Inset: two-dimensional 
Brillouin zone. The F-Y direction corresponds to the [110] 
real-space direction. 

All these results convincingly prove that the Si(l 12)6 X 
1-Ga surface contains a step-edge Ga row and a ter- 
race Ga row. the reduction of the number of partially- 
filled dangling-bonds is the main mechanism determin- 
ing the atomic structure of the metal overlayer. In retro- 
spect, the present findings may have come as no surprise. 
Chemical passivation should always be the key driving 
force for two-dimensional overlayer growth and/or inter- 
face reconstruction. Trivalent Ga is a perfect candidate 
for passivating a silicon surface whereas monovalent met- 
als such as Ag and Au generally induce a more compli- 
cated reconstruction Q,I2>lilEi- This paper, however, 
highlights a more subtle point, that is, the thermody- 
namic and kinetic stability of these and other atom wire 
structures should always be addressed in the context of 
the applicable range of adsorbate's chemical potential. 
Specifically, the fabrication of potentially useful nanowire 
structures on stepped surfaces must take into account na- 
ture's desire to passivate dangling bonds on the terraces. 
The kinetic limitations could be more case-specific. Our 
conclusions should nonetheless have quite general valid- 
ity and therefore reinforces the need for similar analysis 
on other dangling-bond wire arrays, including those high- 
lighted in Ref. |lj , to understand the generic mechanisms 
that control the formation of atomic wires on stepped 
semiconductor surfaces. 



In conclusion, we have analyzed the formation of atom 
wires on vicinal semiconductor surfaces for the case of 
Ga on Si(l 12) . Detailed STM experiments and thorough 
DFT calculations show that the passivation of the sub- 
strate's dangling bonds is the main mechanism determin- 
ing the structure of the overlayer. The Si(112)6 x 1- 
Ga surface consists of two parallel rows of Ga atoms 
adsorbed on the terraces and at the step-edges, inter- 
sected by quasi-periodic vacancy lines. The surface is 
non-metallic and exhibits quasi-one dimensional charac- 
ter in the lowest conduction band. 
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